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Abstract. We investigate solutions for a particular class of linear equations in dendriform algebras. 
Motivations as well as several applications are provided. The latter follow naturally from the intimate 
link between dendriform algebras and Rota-Baxter operators, e.g. the Riemann integral map or Jackson's 
q- integral. 
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1. Introduction 

This paper sets out to develop a systematic study of a particular class of linear equations in Loday's 
dendriform algebras |Lod01| . The guiding principle follows from the intimate connection of such algebras 
to associative Rota-Baxter algebras |AgOO IBax60t IEb02l IRot69j . where analog equations naturally 



appear in the context of applications ranging from numerical analysis to renormalization in perturbative 
quantum field theory. The reader interested in more details may want to consult some of the following 
references [MM IEMPn7b[ IE(;Pn7[ IEMPn7j . 

In the setting of dendriform algebras we obtain the recursive formula for the logarithm of the solutions 
of the two linear dendriform equations Y = 1 + \a ^ Y and Z = 1 — XZ >- a, a € D, in L'[[A]], 
where {D,~<,y) is a unital dendriform algebra. We also present the solutions to these equations as 
an infinite product expansion of exponentials. This way, our work provides a refined approach to the 
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classical Magnus Mag54 and Fer expansions |F58j . well-known in the field of numerical analysis in 



the context of approximations of ordinary first order linear differential and integral equations, e.g. see 
|BCOR98| IIse84( IIse02j . In our approach both expansions involve the pre-Lie and associative products 
naturally associated with the underlying dendriform structure and display their interesting interplay. 
Let us emphasize that the language of dendriform and pre-Lie algebras unveils a hitherto unobserved 
new structure in these expansions, related to operadic aspects in the context of free pre-Lie algebra 
[ChaLivOT] . Moreover, the two linear equations above may be interpreted as first order cases with respect 
to the parameter A. We will introduce equations of arbitrary order in A and show by transforming the 
higher order equations into dendriform matrix equations that its solutions follow from the first order 
ones. 

Putting the two dendriform products together into a single equation: 

(1) X = a + XX yb + Xc^X, 

we find a natural link to Lie bracket flow equations. However, for general dendriform algebras the Lie 
bracket is replaced by a pre-Lie poduct. We will show how the solutions Y = Y(c) and Z = Z{b) of the 
first two order recursions above naturally lead to the solution of the general equation ([1]): 

(2) X = Y * (y-^ y a -< Z~^) * Z. 

Here, the product * is the associative product in the dendriform algebra D. In fact, we observe that in 
the special case where —b = c, the solution X can be represented in terms of two non-commuting group 
actions. 

Eventually our approach brings together the works of Magnus |Mag54 , Fer |F58j and Baxter [BaxGOj 



on linear initial value problems and the corresponding integral equations. A simple example might help 
to elucidate this point of view. Details will be provided in the sequel. Let k he a field of characteristic 
zero and J- an -associative- /c-algebra, say, of operator-valued functions on the real line, e.g. smooth 
n X n matrix-valued functions. Recall that the solution to the initial value problem: 

Y{t) = [Y{t),A{t)] = {YA){t) - {AY){t), t > 0, Y{0) = Yo, 

for A & J^, can be written in terms of the solution of the initial value problem U{t) = A{t)U{t), 
U{0) = 1, i.e. Y{t) = -Adu(t){Yo) =■ -U{t)YoU~^{t). Following Magnus |Mag54| , see also [IMNZnnj . 
we can write U{t) = exp{Q{A){t)), where i}{A){0) = and: 

mm = A Am, 

exp(adn) - 1 

which implies: 

Y{t) = -^dcxp(n{A)(t))(>o)- 
Here, we may remark that Magnus' exponential solution gives rise to the following nontrivial identity: 

Texp(/(A)(t)) := 1 + I{A){t) + I{AI{A)){t) + I{AI{AI{A))) + ■■■ = exp(0(A)(t)), 

I{X){t) = Jq X{s)ds denotes the Riemann integral on T. The expression on the left hand side of the 
first equality is known as the time-ordered exponential, i.e. the T-operation denotes the time ordering 
operator [ORnOj . We may rewrite the Lie bracket fiow equation for A £ J^, as follows: 

(3) Y{t) = [I{Y){t) + Yo,Ait)] = [Yo,Ait)] + ^(t)] = A^it) + 

The crucial observation is that recursion ([3]) actually takes place inside a pre-Lie algebra. Indeed, the 
integration by parts rule implies that the product defined in terms of the Riemann integral and the 
ordinary Lie bracket, i.e. X \> Z := [I{X), Z] satisfies the -left- pre-Lie identity, leading to the pre-Lie 
flow equation: 

(4) Y{t) = Ao{t) + {Y> A){t), Y{0) = Ao{0) = [Yo,A{0)]. 
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One may interpret this as a linear initial value problem on a pre-Lie algebra. When expanding the 
Lie bracket, = I{X)Z — ZI{X), the newly introduced pre-Lie product, X [> Z = [I{X),Z] 

on J^, can be written as a linear combination of two binary compositions, X y Z := I{X)Z, and 
Z ~< X := ZI{X), which we call left- and right-dendriform products, respectively: 

(5) Y{t) = Aoit) + iYyA){t)-iA^Y){t), ¥(0) = Ao{0). 

The rules these left- and right-dendriform compositions have to satisfy so as to provide a pre-Lie 
product are called dendriform axioms |Lod01| . In the next section we will present them in detail. 
Iserles in [IseOlj further generalized the above Lie bracket flow by looking at the initial value problem 
Y = YA- BY, t > 0, y(0) = Yo which can be written: 

Y{t) = YA{t)-BY{t) = {I{Y){t) + Yo)A{t)-B{t){I{Y){t) + Yo) 

(6) = Cit) + iIiY)A)it)-{BI{Ym, 

with C{t) := YQA{t) — B{t)YQ. Using the dendriform products, -< and y, this writes elegantly as follows: 

(7) Y = C + Y y A- B ^Y, y(0) = C(0) = ^0^(0) - 5(0)^0- 

As we will see in the next section, the surprising fact is, that the same dendriform rules also give rise to 
an associative product, and both together, the associative and the pre-Lie product may be used to find 
an exponential solution for Iserles' problem ([7]), and hence also for the above recursions ([1]), respectively 
(0), though in the more general context of dendriform algebras. 

The link to Baxter's work |Bax60j follows naturally from the fact that we solve the above recursions 
using purely dendriform algebra, which allows us to replace the Riemann integral map I by more general 
integral type operators, e.g. Rota-Baxter operators of arbitrary weight |Eb02| . including Jackson's q- 
integral as well as Riemann sums. 

Finally, we will present several applications at the end of the article, including the noncommutative 
generalization of a Spitzer type identity first presented in a commutative context by Vogel in |V63j . 
coming from probability theory. A link to Riccati's differential equation is outlined, which we hope to 
deepen in the near future. We will use planar rooted trees to encode the pre-Lie Magnus expansion and 
also comment on an observation made in an earlier publication |EMj . indicating a reduced number of 
terms in the pre-Lie Magnus and Fer expansions due to the pre-Lie relation. 

The paper is organized as follows. In section [2] we recall the definition of a dendriform algebra and 
introduce the associated left and right pre-Lie products. The augmentation D of a dendriform algebra 
Dhy a unit 1 is also recalled ( |Cha02| . |RonOO| . |Ron02j ). In section [3] we introduce our most general 
equation in D[[A]]: it is said to be of degree (m,n) as it involves monomials of degree < m + 1 (resp. 
n + 1) with respect to the right (resp. left) dendriform product >- (resp. -<). We prove in detail 
that ([2]) gives a complete solution of the case {m,n) = (1,1), which corresponds to the equation ([1]) 
for a,b,c G D, in terms of the solutions of the cases of degree (1,0) and (0,1). The particular (and 
simpler) case, where a = 1 and 5, c G D is also considered. Sections [H and [5] are devoted to the pre-Lie 
Magnus expansion and the pre-Lie Fer expansion respectively, recalling the results of [EMj . In section 
[6] we introduce matrix dendriform algebras, which allow us to solve the equations of degree (2,0) and 
(0,2). The procedure is generalized to equations of degree (m, 0) and (0,n) in the following section. 
Finally, section [8] is devoted to applications in the context of Rota-Baxter algebras, which furnish a large 
collection of dendriform algebras |Eb02] . We touch initial value problems in matrix algebras [IseOlj . an 
identity discovered by W. Vogel [V63] . and show how the Riccati differential equation is linked with 
the degree (2, 0) dendriform equation. Eventually, we introduce a noncommutative Butcher series like 
formula for the pre-Lie Magnus expansion and indicate how the pre-Lie structure may be used to reduce 
the number of terms in it. 
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2. Dendriform power sums expansions 

Loday's dendriform algebra may be seen at the same time as an associative as well as a pre-Lie al- 
gebra. Indeed, it was found that in any dendriform algebra both products can be written as particular 
linear combinations of the dendriform left and right binary compositions. Main examples of dendri- 
form algebras are provided by the shuffle and quasi-shuffle algebra as well as associative Rota-Baxter 
algebras. The link between dendriform algebras and the latter is quite natural as shuffle type products 
and their noncommutative generalizations appear in the construction of free associative Rota-Baxter 
algebras |EG05j . 

We briefly introduce the setting of dendriform algebra. Let A: be a field of characteristic zero. Recall 
that a dendriform algebra |Lod01] over the field /c is a /c-vector space A endowed with two bilinear 
operations, denoted -< and >- and called right and left products, respectively, subject to the three 
axioms below: 

(8) {a^b)-<c = a^{b-<c + b^c) 

(9) {a y b) ^ c = a y {b c) 

(10) ay {by c) = {a ^ b + a y b) y c. 
One readily verifies that these relations yield associativity for the product 

(11) a *b := a ^ b + a y b. 

However, at the same time the dendriform relations imply that the bilinear products > and < defined 
by: 

(12) a \> b := a y b — b ~< a, a <} b := a ~< b — b y a, 
are left pre-Lie and right pre-Lie, respectively, which means that we have: 

(13) (a > 6) [> c - a [> (6 > c) = {b > a) \> c - b o {a > c), 

(14) {a <\b) < c - a < {b <i c) = {a <\ c) < b - a <i {c <i b). 

Observe that a\> b = —b<]a. The associative operation * and the pre-Lie operations >, < all define 
the same Lie bracket: 

(15) l(i,bj ■.= a*b — b*a = a\>b — b\>a = a<\b — b<\a. 
Observe that: 

a*b-\-b\>a = ayb + bya. 

We stress here that in the commutative case (commutative dendriform algebras are also named Zinbiel 
algebras jLod95] . |Lod01| ). the left and right operations are further required to identify, so that a y 
b = b < a. In this case both pre-Lie products vanish. A natural example of a commutative dendriform 
algebra is given by the shuffle algebra in terms of half-shuffles jSch58| . Any associative algebra A 
equipped with a linear integral like map I : A ^ A satisfying the integration by parts rule also gives a 
dendriform algebra, when a ~<b := al{b) and a y b := I[a)b. 

Loday and Ronco introduced in |LR04| the notion of tridendriform algebra T equipped with three 
binary operations, <, > and •, satisfying seven dendriform type axioms: 

(16) {x < y) < z = X < {y -k z), {x > y) < z = x > {y < z), {x -k y) > z = x > {y > z), 
{x > y) • z = X > {y • z), {x < y) • z = x • {y > z), (x • y) < z = x • (y < z), (x • y) • z = x • (y • z), 

yielding associativity for the product x * y := x < y + x > y + x • y. First, observe that the 
category of dendriform algebras can be identified with the subcategory of objects in the category of 
tridendriform algebras with • = 0. Moreover, one readily verifies for a tridendriform algebra (T, <,>,•) 
that {Dt, -<m-,y)i where ■<,:= < + • and ;^:=>, is a dendriform algebra |Eb02j . Any associative algebra 
A closed under a finite Riemann sum operator S : A ^ A satisfying a corresponding modified integration 
by parts rule taking the diagonal into account gives a tridendriform algebra, when a ~< b := aS{b), 
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a y b := S{a)b and a»b := =pa6. Here, the sign of the last composition depends on whether the diagonal 
terms must be subtracted or added, respectively. 

Let A = A(B k.l be our dendriform algebra augmented by a unit 1, such that: 

(17) a ^ 1 := a =: 1 y a l-<a:=0=:a^l, 

implying a*l = l*a = a. Note that 1*1 = 1, but that 1^1 and 1 >- 1 are not defined, see [RonOOj . 
|Cha02] for more details. 

We recursively define the following set of so-called left- respectively right-dendriform iterated elements 
of degree n € N in A[[X\] for fixed xi, . . . , Xn € A: 

w'^\xi,. . . ,Xn) := {■ ■ ■ (Xi y X2) y X3- ■ ■) y Xn 
w'-^\xi,...,Xn) := Xi ^ {X2 < ■ ■ ■ {Xn-1 Xn) ■ ■ ■). 

Defining w^\xi, . . . , Xn) '■= 1 =: Wy\xi, . . . , Xn), we may write these left- respectively right-dendriform 
iterated elements more compactly: 

W^y\xi, . . . ,Xn) := {w^y~^\xi,. . . ,Xn-l)) y Xn 

w^]^\xi,...,Xn) := xi ~< {w':^~^\x2,■■■,Xn))■ 
In case that xi = ■■■ = x„ = a; we simply write w^^^x , . . . , x) = 'w^]^\{x}n) = 'w^^\x) and 

W^^'{X,...,X) = W^^'{{x}n) = W'y'{x). 

Let us recall from Chapoton and Ronco [Cha02l iRonOOl rRon02j that, on the free dendriform algebra 
on one generator a, augmented by a unit element, there is a Hopf algebra structure with respect to the 
associative product ([TT]) . The elements w^'^ := w^\a) generate a cocommutative graded connected 
Hopf subalgebra {H, *) with coproduct: 

(18) A{Wy')=Wy'(E>l + l(!^Wy'+ > ' ^ Wy 

0<m<n 

and antipode S{'w^^) = {—l)^w^^\ It is easy to check that the generate a free associative 
subalgebra of the free dendriform algebra on a for the associative product, so that one can use the 

(n) 

previous formula for the coproduct action on Wy as a definition of the Hopf algebra structure on H. 
As an important consequence, it follows that H is isomorphic, as a Hopf algebra, to the Hopf algebra 
of noncommutative symmetric functions |GKLLRT95j . 

We also define the following set of iterated left and right pre-Lie products (|12p. For n > 0, let 

, . . . , Xyi G A. 

(19) £^''\xi,...,Xn) ■■= {■■■ixi>X2)>X3---)\>Xn 

(20) A"\xi,...,Xn) ■■= Xi < {X2 <\ ■ ■ ■ {Xn-1 <Xn) ■ ■ ■). 

More compactly, for a fixed single element x (z A we may write for n > 0: 

(21) £("+i)(x) = (£(")(x)) > X and r("+^)(x) = x < (r(")(a;)) 

where £^^\x) := x =: A^\x). 

Let us define the exponential and logarithm map in terms of the associative product (jlip . For any 
X G AA[[A]]: 

exp*(x) := 'y^x*"/n! resp. log*(l + x) := — 'y^(— l)"'x*"/n. 

n>0 n>0 
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3. Power sums expansions in unital dendriform algebras 

Let {A, ~<,y) be a dendriform algebra, A its augmentation by a unit 1. Let us introduce the following 
equation of degree {m,n) in A[[A]]: 

m n 

(22) x = a + Y, A^c^55+'^(X, . . . , h,,) + J2 A^a;J+'^(cpi, . . . , Cpp, X), 

q=l p=l 

with ?n, n G N+ and a, bij, I < j < i < m, Cij, l<j<i<nmA. Here, the degree deg{X) € 
of Z is defined as the maximal degree of the left- respectively right-dendriform iterated elements in X 
minus one. In this work we give a detailed account on the solution of the cases (m, 0), (0, n) and (1, 1), 
leaving the general case for further studies. 

We will first explore the equation of degree (1, 1): 

(23) X = a + \J^\x,b) + xJ^\c,X). 
in ^[[A]] for fixed elements a,b,c (z A. 

3.1. Equation of degree (1,1). We first study the case where deg{Z) = (1,1) in (p2|) . Its solution 
follows from the solutions of the cases (1,0) and (0, 1). Later, we will see how the more general cases 
(m,0) and (0,n), m and n bigger than one reduce to the cases (1,0) and (0,1) by embedding the 
dendriform algebra {A, into the matrix dendriform algebra M.n{^) of size N = and N = N^, 

respectively, with 

(m — l)m , (n — l)n 

A^^ := 1 + ^ N^:=l+ ^ ^ ■ 

Relation (j22|) for m = n = 1 simplifies to the following equation: 

(24) X = a + XX yb + Xc^X, 

in A[[X]] for fixed elements a,b,c £ A. 

First, observe that for c = — 6 € ^, relation (p^ reduces to a degree one recursion involving the 
dendriform left pre-Lie product (fT2]) : 

(25) X = a + XX>b. 

in A[[A]] for fixed elements a,b £ A. For this pre-Lie type recursion the special case a = 1 leads to the 
simple solution X = 1, since l[>x = x<\l = 0. 
Let us start to find solutions of ^ 



Proposition 1. Let (A, -<, y) be a unital dendriform algebra. Let Y = Y(c) and Z = Z{b), c,b £ A be 
solutions of the recursions: 

(26) Y = l + Xc'<Y resp. Z = 1 + XZ y b, 

of degree deg{Y) = (0,1) respectively deg{Z) = (1,0), in A[[X]]. One verifies that the product Y * Z 
gives a solution to equation (24^ for the particular case a = 1 £ A. 

Proof. Indeed, a simple calculation shows that: 

Y*Z = {l + Xc^Y)*{l + Zy Xb) 

= 1 + Xc^Y + Z y Xb+{Xc-<Y)*{Z y Xb) 

= l + Xc^Y + Zy Xb + X'^{c^Y) ^ {Z yb) +X^{c-<Y) y {Z y b) 
= 1 + {Xc^Y) {1 + Z y Xb) + {1 + Xc^Y) y {Z y Xb) 
= 1 + {Xc ^Y) -< Z + Y y {Z y Xb) 
= l + Xc^{Y*Z) + {Y*Z)y Xb. 

□ 
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It is obvious that formal solutions to (j26p are given by: 

(27) Y = Y,^'"w'^^\c) resp. Z = ^ A"u;^'^ (6). 

n>0 n>0 

Notice that, due to the definition of the Hopf algebra structure on H, these two series behave as 
group-like elements with respect to the coproduct A, see (llSp . (up to the extension of the scalars from 
k to A; [A] and the natural extension of the Hopf algebra structure on H = ®„>o -f^n to its completion 

H = nn>o with respect to the grading). Hence, we remark here that with respect to equations ([26|) 

the solutions Y = Y{c) and Z = Z{h) of the equations 

(28) Y = 1-\Y c and Z = 1-A6-<Z 
satisfy: 

Y *Y = 1 = Y *Y resp. Z * Z = 1 = Z * Z. 
Indeed we have, using (p7|) and the antipode: 

Y = Y.{-ir\-wt^\c) 

n>0 

= s(5^A"zz;J)(c)) = W = y-i, 

n>0 

and similarly for Z. Hence, we may write Y = and Z = Z~^ , respectively. 

We now slightly generalize equations (|26p by allowing them to start with a non-unital element: 

(29) E = a + \h^E resp. F = c + AF >- d, 

in ^[[A]] for fixed elements a, 6, c, d € A. 

Recall the recursively defined set of left- respectively right-dendriform iterated elements of degree 
n G N in :4[[A]], n > 0: 

w''^'^^\{y]n,x) := y ~< {w':^\{y}n^i,x)), 
w^^^^\x,{y}n) := {w^^~^\x,{y}n^i)) y y, 
in A and for fixed x,y & A. In A[[X]] we find immediately the formal solutions for the equations (f29l) : 

F = a + J^A"-iTz;JHWn-i,a) and F = c + >^''~'wP i^, {d}n-i). 

n>2 n>2 

Our first result is the following: 

Lemma 2. Let A be a unital dendriform algebra. Let Y = Y{b) and Z = Z{d), b,d € A, be solutions 
in A[[X\] to the equations l2^) . respectively. For a,c £ A the equations (2^) are solved by: 

(30) E = Y* (y-^ >- a) and F = (c -< Z'^) * Z 
in ^[[A]], respectively. 

Proof. Let us verify this for the first equation. The other case works analogously. 
Y * (y~i ya) = {1 + Xb^Y)* (y~i >- a) 

= Y-^ ya+ {Xb -< Y) * (y-^ >- a) 

= Y-^ ya+ {Xb ^Y) y {Y'^ y a) + {Xb -< y) ^ (y-^ y a) 
= Y-^ y a + (^{Xb ^ Y) * Y'^^ ya + Xb^(Y* {Y~^ y a)^ 

= y"V a + ((y - 1) * y~^) ^ a + A6 -< (^y * (y^^ y a)) 
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= a + Xb -< (y * (Y~^ ^ a)). 



□ 

An important remark concerning the restriction a, c € A is in order. The equations (|29p reduce to 
the ones in (|26p for a = c = 1. However, for the solutions in (j30p we must be careful with respect to the 
particular properties of the dendriform unit, i.e. we have to avoid 1^1 respectively 1 y 1, see (I17p . 
Indeed, for a G A one may rewrite (I30p : 

E = Y* (y^i ya)=Y*{l^a- A(y~^ y b) y a) = Y * a - XY * ^ b) ^ a. 

Written in this form we may now put a = 1, keeping in mind that b ^ 1. Hence, we find E = Y, i.e. 
that E, with b instead of c, is a solution of the first relation in (j26p as expected. From this we conclude 
our next result, i.e. the solution to equation ([24 



Theorem 3. Let A be a dendriform algebra augmented by a unit 1 (H). Let Y = Y{c) and Z = Z{b), 
c,b & A, be solutions in A[[X]] to the equations Y = 1 + Xc -<Y , Z = 1 + XZ >- b, of degree (0, 1), (1, 0) 
respectively. Then the equation: 

X = a + XX y b + Xc^ X, 
of degree deg{X) = (1, 1) and a ^ A is solved by: 
(31) X = Y * [Y-^ a ^ Z-^) * Z. 

Proof. Let us verify this in detail. Recall equations psp : 
Y * (y~^ ^ a ^ Z"^) * Z = (^(1 + Ac ^ y) * (y~^ >- a< * Z 

= (y-^ ya< z~^) * z + {^{Xc <Y) < (y-i ya^ Z~^)^ * Z 

+( (Xc ^ y) ^ (y~^ >- a ^ z~^)] * z 



(y-^ ya^ Z-^) *Z+\Xc<{Y * {Y-^ ya^ Z"^))) * ^ 

+ (((Ac ^ y) *Y~^) ^ (a ^ Z~^)) *Z 

(y-^ ^ a ^ z-i) * z + (^Ac ^ (y * (y-1 ^ a ^ ^"^))) ^ z 

+(^Ac ^ (y * (y-1 a ^ ^"^))) ^ z + (^((Ac ^ y) * y-^) ^ (a ^ z-^)) * z 

(y~^ ^ a ^ * z + (^Ac ^ (y * (y-^ ^ a ^ >- z 

+Ac ^ (y * (y-^ ^ a ^ z-^) * z) + (((y - 1) * y-^) ^ (a ^ z-^)) * z 

(Ac -< (y * (y-^ ^ a ^ ^ Z + Xc^ (y* {Y~^ ya^ Z~^) * z) + (a ^ Z~^) * Z 

a + Ac ^ (^y * (y~^ ^ a ^ Z-^) * z) + (^a ^ Z"^ + (Ac ^Y) ^ {Y~^ ya^ Z^^)) >- Z 
a + Xc^ (y* {Y~^ ya^ Z"^) * z) 

+ (a ^ Z-i + (Ac ^ y) ^ (y-^ a ^ Z^^)) ^ (AZ y b) 
a + Xc^{Y* (y-^ ^ a ^ Z~^) * z) 

+A(^(a ^ Z"^) * Z + {^[Xc<Y) ^ (y-i ^ a ^ Z"^)) * z) ^ 6 
a + Ac ^ (y * (y-^ ^ a ^ Z"^) * z) + A(^(a ^ Z"^) * Z + (^y ^ (y^^ a ^ Z"^)) * z) I 
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a + Xc^ (y* y Z~^) * + A(^(a -< Z~^) * Z + (y * {Y~^ y Z~^) ] * Z] y b 

-A(((y * Y~^) (a ^ Z-^)) *Z^ yb 
a + Xc^ (y* ^ a ^ Z~^) * z) + x(y * {Y~^ ^ a ^ Z~^) * z) y b 



= a + Xc^X + XXyb 

□ 

Remark 4. Observe that for b,c S A, when using {x y y) ~< z = x y (y ~< z), which imphes 
1 y a ^ 1 = a for a € A, we find: 

X = Y * {Y-^ y a ^ Z~^) * Z 

= y * (a - A(y~^ y c) y a- Xa ^ {b -< Z~^) + A^(y"^ y c) y a < {b < Z~^)) * Z. 

This last expression reduces to X = y * Z, the solution of (j24p . for a = 1. It is tempting in view of this 
observation to set 1 ;^ 1 -< 1 = 1, but it would easily lead to some contradictions. Indeed, if we were 
able to define l^lorl;^lina way compatible with the dendriform axioms, an obvious computation 
leads to 1 )^ 1 ^ 1 = 0! It is then safer to get rid of all expressions involving the undefined 1 ^ 1 or 
1 ;^ 1 in the calculations. 

Corollary 5. Let A be a unital dendriform algebra and a,b € A. Let Z = Z{b) be a solution to 
Z = 1 + XZ y b. The solution to the "left" pre- Lie equation in A: 

(32) X = a + XX>b, 
is given by: 

(33) X = Z-^ *{Zy a-i Z~^) * Z. 

Recall that the solution Z{b) of the equation of degree (1, 0), Z = 1 + XZ ;^ 6 is a group-like element 
with respect to the coproduct A, see (fTSj) . We may therefore write the above solution (j33j) as: 

X = Ad*z-i{Z ya-< Z-^), 

where as usual Ad*p{G) := F * G * F~^. Due to the second dendriform axiom ([9]) we may interpret the 
term Z y a Z~^ as a second group action, which we denote by: 

(34) @z{a) := Zya^ Z~\ 

Lemma 6. Let A be a unital dendriform algebra and a,b G A. Let A, B be two solutions of the equation 
of degree (1,0). 

(1) GaoGb(x) = eA*B(x) 

(2) Qa o e^-i(x) = e^-i o 6^(2;) = X. 

Proof. The second property follows from the first. Let us verify the first identity: 

Qa°Qb{x) = Ay {B y X < B-^) < A-^ 

= {A*B)y X <{A*By^ = Qa*b{x). 
This follows immediately from the dendriform axioms. □ 

We remark that the two group actions in general do not commute, i.e. Ad*j^ o 0^ ^ o Ad*j^. 

Recall that the Dynkin operator is the linear endomorphism of the tensor algebra T{X) over an 
alphabet X = {xi, . . . , . . .} into itself the action of which on words yi ■ ■ .yn, 2/i € X is given by the 
left-to-right iteration of the associated Lie bracket: 

D{yi, . . . ,yn) = [• • • [[yi,y2],y3] ■ ■ ■,yn], 
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where [x, y\ := xy — yx. The Dynkin operator is a quasi-idempotent: its action on a homogeneous 
element of degree n satisfies = nD. The associated projector D/n sends Tn{X), the component of 
degree n of the tensor algebra, to the component of degree n of the free Lie algebra over X. The tensor 
algebra is a graded connected cocommutative Hopf algebra, and it is natural to extend the definition of 
D to any such Hopf algebra as the convolution product of the antipode S with the grading operator A^: 
D := S-kN |PK02l IEGP061 IEGP071 IEMP07j . This applies in particular in the dendriform context to 
the Hopf algebra H introduced above. We will write Dn for Dop^, where pn is the canonical projection 
from T{X) (resp. H) to Tn{X) (resp. if„). 

Lemma 7. |EMP07j For any integer n > 1 and for any x € A we have: 



Let us come back to equation (p^ respectively CorollaryO Then, using D{1) = 0, we see immediately 
that the solution Z = Z{a) to the equation Z = 1 + \Z ^ a is mapped to the solution 

X := X^^D{Z) = Ad*z-i oQz{a) = Z~^ * {Z y a -< Z'^) * Z 

of the corresponding pre-Lie equation (j25p for a = b = —c € A: 

X = a + XX \> a. 



In this section we recall results from |EMj where we obtained a recursive formula for the logarithm 
of the solutions of the two dendriform equations X = 1 + Xa ^ X and Y = 1 — XY ;^ a in ^[[A]], a e A, 
of degree (0, 1), (1,0), respectively. 

Let us introduce the following operators in (^4, ~<, y), where a is any element of A: 

L^[a]{b) := a -< b Ly[a]{b) := a y b R^[a]{b) := b -< a R^[a]{b) := b y a 
L^[a]{b) := a<\b Lf^[a]{b) := a \> b R^[a]{b) := b <i a Ri^[a]{b) := b > a 

In the following theorem we find a recursive formula for the logarithm of the solutions of the equations 
in ()26p . We will call it the pre-Lie Magnus expansion for reasons which will become clear in Remark [9l 

Theorem 8. |EMj (pre-Lie Magnus expansion) Let Q,' := ^'{Xa), a ^ A, be the element of XA[[X]] 
such that Y = exp*(0') and Y~^ = exp*(— $7'), where Y and Y~^ are the solutions of the two equations 
Y = \ + Xa ^Y and Y~^ = 1 — XY~^ y a of degree (0, 1), (1,0), respectively. This element obeys the 
following recursive equation: 



(35) 



D{w'^\x)) = 




4. The pre-Lie Magnus expansion 



(36) 




or alternatively: 



(37) 




where the Bi 's are the Bernoulli numbers. 



Returning to Theorem [3] we see that the solution of equation ()24p writes: 

X = exp*{n'{c)) * (exp*(-!7'(c)) ^ a ^ exp* (0'(-6))) * exp* (-0'(-6)) 
and the solution to the pre-Lie recursion, i.e. equation (j24p for c = —b: 

^ = ^'^cxp*(~n'(-b)) ° 0cxp*(f7'(-fe))(a)- 
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Remark 9. (classical Magnus expansion) Let T be any algebra of operator-valued functions on the real 
line, e.g. smooth n x re matrix- valued functions, closed under integrals I{U){x) := /g U{s)ds. Then, 
is a dendriform algebra for the operations: 

X X 

A -< B{x) := A{x) ■ j B{t)dt A y B{x) := j A{t)dt ■ B{x) 



with A,B^T. This is a particular example of a dendriform structure arising from a Rota-Baxter 
algebra structure -of weight zero |AgOO| . We refer to Section [8] for further details on Rota-Baxter 
algebras and their connections to dendriform algebras. Here, let us simply mention that the Rota-Baxter 
operator on giving rise to the dendriform structure is the Riemann integral map: I{A){x) := A{t)dt. 
In this particular example we see that the pre-Lie Magnus expansion (j36p with the pre-Lie product 
L^[A]iB) = {A> B){t) = [I{A),B]{t): 

(38) Q'{A){s) = n{s) = ^^d'-^llM^s)), 

m>0 



coincides with Magnus' logarithm Mag54| , ^}{t) = \og{Y(t)), of the solution of the initial value problem 



Y{t) = A{t)Y{t), y(0) = Yo = l, such that: 

y(t) =exp( /"*0(s)ds) = 17(0) = 0, 

■^0 n>0 

for more details on Magnus' work see also |Cals06l IGKLLRT951 UNMl IKU891 iMPTOl [5t87l IWilGTl [Z96] . 
By now we hope the above used prime notation, Q', in equation (p6]) and ([37|) of Theorem [8] has become 
self-evident. Notice that if we suppose the algebra J-" to be unital, the unit (which we denote by 1) 
has nothing to do with the added unit 1 of the underlying dendriform algebra. In fact, we extend the 
Rota-Baxter algebra structure to by setting: 

(39) R{1) := 1, ^(1) := -1 and 1.x = x.l = for any x eT. 

This is consistent with the axioms (jl7p which in particular yield 1 y x = R{l)x and x ^ 1 = —xR{l). 
Later, in Section [8] we will provide more details in the context of applications. 

Returning to the two group actions, we describe the infinitesimal actions. Let ^4 be a unital dendriform 
algebra. Let Y = exp*(r2'(Aa)), a £ A, he the solution of the equation Y = 1 + Xa ~<Y oi degree (0, 1) 
mA[[X]]. 

(1) @y{x) =Y>- X '<Y-'^ r^x + \a>x + 0(A2) 

(2) Ad*Y{x) = y * x * ~ X + X[a, x] + 0(A2). 

5. The pre-Lie Fer expansion 

Following Fer's original work |F58j . see also |Ise841 lIMNZOOt IWil67| . we now make a simple Ansatz 
for the solution of the degree (0,1) equation Y = 1 + Xa ~<Y : 

Y = exp*{Xa)*Vi 

and return this into the recursion. This then leads to what we named the pre-Lie Fer expansion. 

Theorem 10. |EMj (pre-Lie Fer expansion) Let {A, ~<, y) be a dendriform algebra augmented by a unit 
1 (7?). Let Uq := Xa, := U^{a), n G N, a G A, be elements in XA[[X]], such that 



n>0 n>0 
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where Y and Y~^ are solutions ofY = l + Xa^Y and Y~^ = 1 — XY~^ >- a, respectively. Then these 
elements obey the following recursive equation: 

-lyi 



(1 + 

The presentation given here in the context of dendriform algebras, using the corresponding pre-Lie 
products, reduces to Fer's classical expansion |F58j . see also |MZ97j . when working in the dendriform 
algebra introduced in Remark [H 

6. An equation of degree (2,0) 

Let D he a dendriform algebra and D its augmentation by a unit 1, see ()17p . We are interested in 
solving the following equation of degree (2,0) in -D[[A]]: 

(40) X = d + XXyc + X'^{Xyb)ya, 



for a,b,c € D and d ^ D. We introduce matrix dendriform algebras, in which equation (j40p will be 
transformed into a system of order (1,0) equations. 

6.1. Matrix dendriform algebras. Let (Z), -<, >-) be a dendriform algebra. Then the space M.n{D) 
of square n x n-matrices with entries in is a dendriform algebra, with operations defined by: 

n n 

(a ^ b)ij := ^fli/fc ^ bkj, (a >- b)ij := ^ajfc >- bkj. 

1=1 i=l 

We can indeed verify the first dendriform axiom: 

n 

[{a ^b) ^ c]ij = ^(a -< b)ik -< Ckj 

k=l 

n 

k,l=l 

n n 

= au ^ {kk * Ckj) = '^aii ^ {b* c)ij = [a -< (6 * c)]ij, 

k,l=l 1=1 

hence {a ~< b) ~< c = a ~< {b * c) in A^„(D), and the two other axioms can be checked similarly. If 
D := DQ)k.l (resp. Mn{D) := A4„(L') © A;.l„) is the dendriform algebra (resp. 7W„(Z))) augmented 
with a unit, then we can identify the unit 1„ with the diagonal matrix with dendriform units I's on the 
diagonal and elsewhere. 

6.2. Transformation into a dendriform system of degree (1,0). First, we will show how to 
transform equation (|40p for d = 1 into a "first-order" equation in A^2(-C)[[A]]. Observe the size of the 
matrices, i.e. 2 = 1 + (2 - 1)2/2. 



Lemma 11. Let Z G Al2(-D)[[A]] be the solution of the "matrix dendriform" equation of degree (1,0).- 

(41) Z = l2 + XZ y M. 

Then for any v = {vi, V2) G k"^ the solution in (D[[A]])^ of the equation: 

(42) Y = vl2 + XY y M, 

with Y y M := (Yiy Mu +Y2>~ M21, Yi y Mu + I2 ^ M22) is given by: 

vZ := {viZii + V2Z21, V1Z12 + ^2-^22)- 
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Proof: We immediately check: 

vZ = VI2 + Xv{ZyM) 
= VI2 + X{vZ) y M. 

□ 

Proposition 12. Let X be the solution of equation in D[[X\] with d = 1. Then the vector Y := 
{X, XX >- b) in (Z)[[A]])^ is the solution of equation with v = (1,0) € k'^ and: 

^=(" n 
\a 

Proof: One easily checks that: 

VI2 + XY y M = (1, 0) + X{X, XX yb)y (^^ 

= {l,0) + X{Xyc + X{Xyb)ya,Xyb) 
= (1, 0) + {X - 1, XX y b) 
= Y. 

□ 

Using Theorem [HI i.e. the pre-Lie Magnus expansion, we find immediately the solution to (|4ip : 

Corollary 13. Let il.' = log* Z S A^2(^[['^]]); where Z is the solution of equation |^ j[ j with M 
'c 6^ 



^ Q . , given by the pre-Lie Magnus expansion: 

n'{XM) =XM - ^X^M [> M + ^X^{M t> M) >M + X^^M t> {M \> M) ■ ■ ■ . 

Then the solution X of equation d = \, is such that the line vector {X, XX y b) is the first line 

of the matrix Z = exp*(0'), i.e.: 

X = (1,0)Z(1,0)*. 

The following corollaries are readily verified. 

Corollary 14. Let X be the solution of equation in D[[X]]. Then the vector Y := {X, XX y b) in 
(Z)[[A]])^ is the solution of equation: 

(43) Y = v{d* I2) + XY y M, 

with V = (1,0) € A;^ and M as in Proposition[ 



Here d* I2 is the diagonal matrix in A^2(-C) with d on its diagonal and else. Using Lemma [2] we 
immediately get: 

Corollary 15. Let fl' = log* Z G A42(-D[[A]]); where Z is the solution of equation (f^ip with M as in 
Proposition MSX given by the pre-Lie Magnus expansion: 

n'{XM) =XM - ^A^M > M + ^A^(M o M) \> M + X^^M t> (M t> M) ■ ■ ■ . 
The line vector {X, XX y h) is the first line of the matrix U := [dl2 -< exp*(— 17')) * exp*(r2'). 

The equation of degree (0, 2) in -D[[A]]: 
(44) X = l + Xc^X + X^a-<{b<X). 

can be treated similarly: 
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Proposition 16. Let X he the solution of equation in D\\X\\. Then the column vector Y := 

{X, X b ^ Xy in {D[[X]])^ is a solution of equation: 

(45) Y = l2V + X ^Y, 

with V = (1,0)* G k'^ and: 




Proof. The solution of equation in D[[X]] is exactly the solution of equation ([lO]) in I?'[[A]] where 
D' is the dendriform algebra (D, ^, with x ^ y := y y x and x ^ y := y ~< x. □ 



7. Higher order equations 
The general degree (m, 0) equation writes: 



(46) 



-^^=000 + ^ X''uJ^y^^\X, Ogl,..., aqq). 

q=l 



The cases m = 3 and m = 4 can be worked out by hand and will give a hint for the general case. Indeed, 
the length 4 line vector: 

Ys := {X, XX y 021, XX y a^, X'^{X y asi) y 032) 
is the solution of the equation: 



Y3 = (aoo, 0, 0, 0) + XY3 y M3, with M3 



The length 7 line vector: 

, AX y a2i, XX y 031, 
is the solution of the equation: 



Y4 = (aoo, 0, 0, 0, 0, 0, 0) + XY4 y M4, with M4 



We recognize the matrix as the upper left 4 x 4-submatrix of The solution of the general 
equation (j46p will be the first coefficient of a vector Y^ of length 1 + m{m — l)/2, whose coefficients 
(discarding the first one) are given by X^u;i^\x, aq\, . . . ,aqj), 1 < J < g < m. They are displayed 
according to the lexicographic order of the pairs {q,j) indexing them. The vector Y^ is the solution of 
the equation 

Yn, = (aoo, 0, . . . , ) + XYmy M^, 
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where the square matrix of size 1 + m(m — l)/2 is recursively defined by: 



(47) 



Mm-i Ai 
A2 A3 



/ ami 




m—l 



/ 













0\ 




0/ 



V 



























W 








\ 




'^m(r?i— 1) 





We proceed then as in paragraph 16.21 according to Lemma [2j The line vector Ym is the first line of the 
matrix of size := 1 + m(m — l)/2 defined by ((aoo * Itv) -< Z*~^^ * Z, where Z is the solution of the 
order (1,0) matrix dendriform equation: 

(48) Z = lN + ZyM,n- 

In the simpler case when aoo = 1 instead of being an element of D, the line vector is just the first 
line of the solution Z of equation 



The solution of general equation of order (0, n) can then be derived immediately, by the same trick as 
in the proof of Proposition [T6j Higher-order equations of type (|22p of order (m, n) remain to be solved. 

8. Applications 

For completeness let us recall the recursions and its solutions we described in the foregoing sections. 
Let (D[[A]], -<, >-) be a unital dendriform algebra, a,b £ D. The principal equations are: 

Y = l + Xa^Y Z = 1- XZ y b, 



of degree (0,1) and (1,0), respectively, and with formal solutions Y = Y{a) = ^„ 



>o 



(n) 



(Aa) and 

Z = Z{b) = En>o^>-"^(--^^) ill D[[>]]- RecaU that Z{a) = Y~^{a) gSI- In Theorem [8] (pre-Lie Magnus 
expansion) we showed that Y{a) = exp*{W{Xa)), whereas Theorem [10] (pre-Lie Fer expansion) concludes 

that Y{a) = nn>o ^^P*(^n('^))) ^0 ~ Another pair of recursions in -D[[A]] are: 

E = a + Xb^E F = c + Fy Xd, 

a,b,c,d G D. Lemma [2] provides the solutions to these equations in D[[A]]. They can be expressed in 
terms oiY = Y{b) and Z = Z{d), i.e. E = Y * [Y'^ y Aa) respectively F = (Ac -< Z'^) * Z. Finally, 
putting together the former two recursions we arrive at the equation of degree (1, 1): 

X = a + XX y b + Xc ^ X, 

a,b,c D, and its solution in terms of 1" = Y{c) and Z = Z[b): 

X = Y * (y~^ ^ a -< Z'^) * Z. 

For c = —b £ D we arrive at the pre-Lie recursion in D: 

X = a + X > b, 

and its solution in terms if Z = Z{—b): X = Ad*^_i o 6^(0). 

Recall from the introduction the link between recursions in dendriform algebras and simple initial 
value problems, say, for instance in algebras of matrix valued functions closed under the Riemann 
integral map. In fact, the link here is provided by the relation between associative Rota-Baxter and 
dendriform algebras. 
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8.1. Rota— Baxter algebras. Recall |Bax601 IEb02l IRot69j that an associative Rota-Baxter algebra 
(over a field k) is an associative A;-algebra A endowed with a /c-linear map R : A ^ A subject to the 
following relation: 

(49) R{a)R{b) = R{R{a)b + aR{b) + Oab). 

where 9 (z k. The map R is called a Rota-Baxter operator of weight 9. The map R := —9id — R also 
is a weight 9 Rota-Baxter map. Both the image of R and R form subalgebras in A. Associative Rota- 
Baxter algebras arise in many mathematical contexts, e.g. in integral and finite differences calculus, 
but also in perturbative renormalization in quantum field theory |EGM06] . 

A few examples are in order. First, recall the classical integration by parts rule showing that the 
ordinary Riemann integral is a weight zero Rota-Baxter map. Moreover, on a suitable class of functions, 
we define the following Riemann summation operators: 

[x/e] ix/e]~i 

(50) Re{f){x):=Y,()f{n0) and R'eiDix) := ^ 9f{n9), 

n=l n=l 

which satisfy the weight —9 and the weight 9 Rota-Baxter relation, respectively. The summation 
operator: 

(51) 5(/)(x):= J]0/(x + 0n). 

n>l 

also satisfies relation (149 p on a suitable class of functions. Another useful example is given in term of 
the g-dilatation operator, agf{x) := f{qx). In fact, one shows that the g-summation operator: 

(52) S,{f){x) := 5]</(x) = ^/(<z"x) 

n>0 n>0 

satisfies relation ()49p on a suitable class of functions, with 9 = —1. In fact, one may think of (j49p as a 
generalized integration by parts identity. Indeed, the reader will have no difficulty in checking duality 
of (09]) with the 'skewderivation' rule |CEFG06j : 

d{fg) = d{f)g + fdig) + 9d{f)d{g), 

i.e. Rdx = X + c and dRx = x, x ^ A and c G ker d. For instance, in the case of the g-summation 
operator ([52I1 . we find dqf{x) := {id — aq)f{x) = /(x) — f{qx), satisfying: 

dqifg) = dq{f)g + aq{f)dqg 

= dq{f)g + fdq{g)-dq{f)dq{g). 

Prom this one can prove the Rota-Baxter property for the (/-summation operator Sq (weight 9 = — 1), 
which is equivalent to the so called (^-integration by parts rule: 

S,{f)Sq{g) = Sq{Sq{f)g) + Sq{aq{f)Sq{g)). 

We call a A;-algebra A with both, a Rota-Baxter map R and its corresponding skewderivation d a Rota- 
Baxter pair [A,R,d). Although, we should underline that, in general, Rota-Baxter algebras do not 
come with such a skewderivation, e.g. let A be a K-algebra which decomposes directly into subalgebras 
Ai and A2, then the projection to A\, R : A A, R{ai,a2) = ai, is an idempotent Rota-Baxter 
operator, i.e. of weight 9 = —1, see |CEFG06t lEMj for more details. 

Proposition 17. |Eb02j Any associative Rota-Baxter algebra gives a tridendriform algebra, {Tfi,<,> 
,•0), in the sense that the Rota-Baxter structure yields three binary operations: 

a < b := aR{b), a > b := R{a)b, a»e b := 9ab, 



satisfying the tridendriform algebra axioms (11) 
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The associated associative product *g is given by 

(53) a*eb:= aR{h) + R{a)b + Oab. 

It is sometimes cahed the "double Rota-Baxter product", and verifies: 

(54) R{a*gb) = R{a)Rib), R{a *e b) = -R{a)R{b) 

which is just a reformulation of the Rota-Baxter relation (j49p . In fact, one easily shows that for any 
Rota-Baxter algebra A, the vector space A equipped with the new product (j53p is again a Rota-Baxter 
algebra, denoted Aji, though in general non-unital. 

From the general link between dendriform and tridendriform algebras it simply follows: 

Corollary 18. Any associative Rota-Baxter algebra gives rise to a dendriform algebra structure, {Dr, -< 
, >-), given by: 

(55) a -< 6 := aR{b) + Bab = -aR{b), a y b := R{a)b. 
The dendriform pre-Lie products, (jl3|) and (jl4p . can be written explicitly: 

a>b= [R{a), b] - 9ba a<ib= [a, R{b)] + 9ab. 

Recall |EMj that we extend the dendriform structure {Dr, ~<, y) on the Rota-Baxter algebra {A,R) 
to -Dr by setting R{1) := 1, -R(l) := —1, and 1.x = x.l = for any x G Dr. This is consistent with 
the axioms (jl7p which in particular yield 1 y x = R{l)x and x ^ 1 = —xR{l). 

8.2. A link with generalized initial value problems. Let us go back to Remark [9l Assume for 
the moment that we work in a noncommutative unital function algebra J-', say, suitable n x n-matrix 
valued functions, with a Rota-Baxter pair {R,d). Let us consider the initial value problems (IVPs) 

(56) Y ■.= dY = YB - CY, Y{0) = Yq 

with B,C ^ J- . Iserles in |Ise01j looked at this equation in the classical setting of Remark El where its 
straightforward solution is given in terms of the classical Magnus expansion (j38p : 

(57) Y = -exp(0(C))yoexp(-J7(S)). 

Here U := exp(0(5)) and V := exp(-0(C)) solve the IVPs, U = -UB, U{0) = 1 and F = CV, 
V{0) = 1, respectively. 

Coming back to equation (j56p in the context of a Rota-Baxter pair {R, d) of weight zero, we may 
transform it into the recursion for Y: 

Y = A + R{Y)A- BR{Y), A = YoB-CYo, 

which now becomes a dendriform equation of degree (1, 1) for the element X := Y: 

X = A + X y B + C -< X. 

One verifies by using the relation between d and R as well as (j54p . that R maps the solution X, see 
pTl). of this equation to i^, i.e. R{X) = Y. 

However, in the light of the fact that Rota-Baxter algebras in general may not come with a corre- 
sponding skewderivation, it is crucial to realize that X is simply an element in Ar, the Rota-Baxter 
algebra with product (I53p . Indeed, we would like to emphasize that our findings in the context of 
dendriform algebras hold through for any associative Rota-Baxter algebra {A,R). 

Hence, starting with a Rota-Baxter algebra {A,R), Theorem [8] (resp. Theorem llOp gives us the 
solution Y = exp(i2(r2'(a)) to the equation Y = 1 + XR{aY). We may interpret this as the "integral- 
like" equation corresponding to the IVP, Y = aY, Yq = 1. Though, we only used (I54p as well as 
P7p . and interpret Y purely as an element in An- Theorem [3] then implies for 6 = the solution to 
recursion E = Xa + XR{cE), given by = exp{R{p,' {c))R{exp{R{^l' {c))a) , analogously for c = 0. From 
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a dendriform algebraic point of view the classical IVP (j56|) should be replaced by an equation in terms 
of the pre-Lie product: 

X = a + X [>b, 

a,h ^ A, where its solution (j33p again is an element in A^: 

^ = ^^exp*(-f7'(-6)) ° 0exp*(f7'(-fe))(a). 

The formal limit of ^ reduces to the Lie bracket flow like equation. It is important to remark that, 
as there are two natural Rota-Baxter maps on A, i.e. R and R, we can map X G Aji either via R to 
Y ■= R{X) G A or via RioY := R{X) E A. 

8.3. A link with the Riccati diflferential equation. The Riccati differential equation writes: 

(58) y=p + qy + ry^, 

where y,p, q and r are differentiable functions on the real line. We suppose for convenience that r never 
vanishes. Looking at the transformation: 

(59) y = — 

rw 

one easily checks that the solutions of (I58p are also solutions of the following linear homogeneous 
second-order differential equation: 

(60) riv — (r + qr)w + pr'^w = 0. 

Now consider equation (j40p . Recall Remark [9l in this particular situation, we find a dendriform algebra 
structure on the set J- of twice differentiable functions on M, which we denote by D, with dendriform 
left and right products: 

(61) f^g:=f-R{g) f y g := R{f) ■ g. 
Here R = I is the Riemann integral operator on defined by: 

(62) Rif){x):= r f{t)dt. 

Jo 

Recall from |EMP07bl IEM] that the operator R (which is a Rota-Baxter map of weight zero) extends 
to D by putting -R(l) = 1, where 1 is the usual unit of the algebra of functions, i.e. the function equal 
to 1 everywhere. The operator R extends naturally to i5[[A]] by A-linearity. Starting from a solution X 
of ([I0|1 and putting y := R{X) we have: 

(63) y = l + XR{R{X)c) + X^R(^R{R{X)a)by 
and hence: 

(64) y = \yc + \^R{ya)h. 

We suppose at this point that 5 is a nowhere vanishing function on M. Dividing out by b and differen- 
tiating once again we obtain: 

fn^\ I" ^ ■ , d c c. 2 

(65) 6^-6^^ = ^^d^6 + V + ^^"' 
hence: 

(66) by-{b + Xbc)y-iX~ + X^a)b^y = 0. 

We recognize a formal version of equation (I60p . with p = — (^f + A^o), q = Ac and r = b. 
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8.4. Vogel's identity. We present here a generalization of an operator identity and its solution, found 
by Vogel in |V63| in the context of probability theory, to noncommutative associative Rota-Baxter 
algebras. 

Let {A, R) be an associative Rota-Baxter algebra of weight 6 (^K, a,b,c € A. The following equation 
defined in ^[[A]] appeared in |V63j : 

(67) X = a + XR{X)b + \cR{X). 

For A being commutative Vogel found the solution, X = B — C , where C = C{b, c) and B = B{b, c) are 
defined in ^[[A]] as follows: 



(68) C := exp(^R{a{X))jR[aexp{R{9~Hog{l + Xec))+R{9~^log{l + Xeb)) 

(69) B := exp(^-R{a{X))^R(aexp{R{e~^log{l + X9c))+R{e~Hog{l + Xeb)) 
and 9a{X) := log(l + X9b) — log(l + X9c) 0. It is obvious that Vogel's identity simply translates into the 



dendriform identity (124^ : 

X = a + XX yb- Xc-iX, 
with its general solution X = Y * {Y~^ y Z'^) *Z,Y = Y{-c) and Z = Z{b) = y(-6)~\ i.e.: 

X = exp*(!^'(-Ac)) * (^exp*(-f]'(-Ac)) ^ a ^ exp* ((f^'(-A6)))) * exp* (-f]'(-A6)) . 

Rewriting it into the Rota-Baxter representation this gives the solution to Vogel's identity in a non- 
commutative Rota-Baxter algebra. It is an interesting exercise to retrieve Vogel's solution from the 
dendriform one in the case of a commutative Rota-Baxter algebra. Hence, first observe that in this 
case: 

a^b = a)^b-b^a = R{a)b - {-bR{a)) = R{a)b - 9ba - bR{a) = [R{a), b] - ba = -9ba. 
Such that: 

n'iXa) = y ^L^[n'](-){Xa) = y ^[-9rx{Q'{Xa)ra = 

^-^ nl ^-^ n! 9 

n>0 n>0 

By recalling our convention (j39|) and using commutativity as well as the well known identity for the 
product §^,9a*b = R{a)R{b) - R{a)R{b), we find: 

,~,log(l + A^c),. , ~.log{l + X9b)^.~f ,^.log{l + X9c). ~.log{l + X9b).. 
X= exp(i?(^^-^ '-)) exp(-i2(^^-^ L)^ Ri^aexp{R{^^^ + «(^^-^ -)) 

, p.log(l + A£c) . ^p^Ml + A£6) . / ^^^ log(l + Agc) ^ , -^ log(l + Ag5) . 
-exp(-i?( ))exp(i?( ))R[aexp[R{ ) + R{ ^ 

8.5. Rooted trees and the coefficients in the pre-Lie Magnus expansion. The final section 
is devoted to a closer analysis of the pre-Lie Magnus expansion. Iserles and N0rsett |IN99j used pla- 
nar rooted binary trees to disentangle the combinatorial structure underlying the classical Magnus 
expansion, see also |IMNZOO| IIse02j . Here we propose the use of planar rooted trees to represent the 
pre-Lie Magnus expansion. This, of course, incorporates the classical case but, is more in the line of 
-non-commutative- Butcher series. 

However, let us first use planar rooted binary trees. Recall rt'{Xa) = X^„>o -^"^^nC*^)- 

Bn J (m) rp./i 



This writes out as: 



n>0 



$7'(Aa) = Xa- X^-a \> a + A^^(a \> a) \> a + X^—a > (a > a) - A^-((a \> a) \> a) \> a 
2 4 12 8 



^We should remark that, following our convention, Vogel uses a commutative Rota-Baxter algebra of weight 9 — —1. 
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— -^"^^ (('^ \> {a[> a)) 0> a + a[> {{a[> a) [> a) + {a[> a) [> {aO> a)j + . . . 
At fifth order we have ten terms. Using planar binary rooted trees to encode the pre-Lie product: 

~ a O a, 

we find at fourth order: 

-I A - UA + A + 

In [EMj we observed that from order four upwards the left pre-Lie relation (jl3p implies a reduction in 
the number of terms. Indeed, at order four we observe a reduction from the four terms above to the 
following two terms: 

"eA ~ 12 

thanks to the left pre-Lie relation: 



To go beyond this order we propose another approach in the spirit of non-commutative Butcher series. 
In fact, one can write the pre-Lie Magnus expansion like a classical Butcher series using ordinary planar 
rooted trees. By definition a planar rooted tree t is made out of vertices and nonintersecting oriented 
edges such that all but one vertex have exactly one incoming line. For a tree t we denote by V{t) and 
E{t) the sets of vertices and edges, respectively. By f{v), v G V{t) we denote the number of outgoing 
edges, i.e. the fertility of the vertex v of the rooted tree t. The degree of a tree deg{t) is given by its 
number of vertices. We denote the graded vector space of planar rooted trees by Tpi := ©m>o^;'^^- Let 
t = B^{ti, ■ ■ ■ ,tn), where denotes the grafting operation, joining the n roots via n new edges to a 
new vertex, which becomes the root of the new tree t and deg{t) = deg{ti) + 1. Now, we define 

the linear map a : Tpi ^ k: 

-w-tn-c-)- n ^- 

1=1 v€V{t) ■' ^ ' 

We see immediately that this function maps a rooted tree t containing vertices of fertility f{v) = 2n+ 1, 
n > 0, V G V{t) to zero. 



Lemma 19. ker{a) = {t eTpi\3v e V{t) : f{v) = 2n + l,n > 0} 

^pl 



In the following we denote by 73^ the graded vector subspace of planar rooted trees excluding trees 



t with vertices of fertility f{v) = 2n + 1, n > 0, v £ V{t). Recall r^^\ai, . . . ,an) := ai 1> (a2 > (03 > 
• ■ • (dn-i l> o-n)) ■ ■ ■ ) • ^ext, we introduce the tree functional F, F[. ](a) = a and for t = B^{ti • • ■ 
n > 1: 

F[t]ia) :=r^^^'\F[h]ia),...,F[tn]ia),a). 
Theorem 20. The pre-Lie Magnus expansion can be written: 

(70) n\a)= a{t)F[t]ia). 

pi 

Proof. The fact that the sum only includes planar rooted trees with even vertices and vertices of fertility 
one follows form Lemma [19l We just have to prove that the RHS of (j70|) verifies equation ([37|) (with A 
set to 1). Indeed, if we denote by this right-hand side we have: 

n" = Y '^{t)F[t]{a) 

pi 
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= E E ^c^{ti)---a{trn)r^^^'\F[h]{a),...,F[tn]{a),a) 
= E|j^£"^'^(21;;ifi:'«)' 

"^-•^ m times 

hence n" = n'. □ 

Remark 21. Let T{z) = 'Ylk>o'^kZ^ be the Poincare series ofT^^'. Using the methods of |D04j we can 
prove that T is the solution of the following functional equation: 

(71) T{z) - z^T{zf = r^^, 
which in turn yields the recursive definition of the coefficients T„; 

(72) Tn = l+ ^f^-?^- 

p,q,r>0,p+q+r=n—2 

This sequence matches with sequence [A04-9130] in Sloane's encyclopedy [S]. On can see this by noticing 
that the series W = zT satisfies the functional equatioi^: 

(73) W{z) - W{zf 



1-z' 
z{ z-l){z + l ) 
z^ - z + 1 
1, 1, 2, 4, 10, 26, 73, 211, 630, 1918, .. 



hence W is the reversion of the series z i— > ^ . The first terms of the sequence are 

yO -.1-1 



This simple formula provides a simple and efficient way to expand Magnus' -pre-Lie- recursion. Let 
us check a few examples: 

Sl'(a) = a(. )F[. ]{a) + a(I )F[l ]{a) + a (i)F[i](a) + a(A)F[A](a) + --- 

= a + B^rg\F[.]{a),a)+Blr^^\F[:]{a),a) + ^rg\F[.]{a),F[.]{a),a)+... 

Recall that here Tpi denotes the vector space of planar rooted trees. Hence, at order four we find: 

n'ia) = ••• + «(! )F[\ ]{a) + q(A )F[A ](a) + afi )F[A ](a) + a(A )F[k ](a) + • • • 

= --- + BfA^\F[l]{a),a)+^B,r^^\F[:]{a),F[.]{a),a)+^B,rl^\F[.]{a),F[:]{a),^ 

+ Bi|^rg)(F[A](a),a)+.... 

Recall that B^ = 0. Due to the recursive nature of the pre-Lie Magnus expansion at order five the 
following set of 8 trees appear: 

{B+{\), S+(A), B+{.,. ,.,.), B+{A,.), i?+(.,A), i?+(I,I), S+(i,.), i?+(J)}. 

Of course, the first two terms have different coefficients. Other trees do not appear due to the fact that 
odd Bernoulli numbers are zero, i.e. -B3 = 0. Hence, for the order five contribution we find: 

^5(0) = ((('^ > a) {> a) [> a) > a - [a > ((a > a) > a)) > a 



"^We greatly thank Loi'c Foissy for having indicated this point to us 
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B B"^ 

+ ^ at> [at> {a\> {at> a))) + ^ (a [> (a > a)) [> (a > a) + 

B^ B 
+ ^ a [> ((a > (a t> a)) [> a) + ^-B^ (a t> a) > ((a > a) > a) 

-B 

+ -^BI ((a > a) > a) (> (a [> a) + a > (((a > a) > a) [> a) . 

Using the left pre-Lie relation we may write: 
-^BI ((a > a)> a)\> {a\> a) 

B B B 

= -^B'1 a \> (((a \> a)\> a)> a) - -^Bf [a \> {{a [> a) [> a)) [> a + -^Bf (((a \> a) \> a) \> a) \> 

leading to seven terms: 

B2 5 B2 1 

^5(0) = ~Bi-^-^ (((a > a) > a) [> a) > a - ^1^2 ^'^ ^ '^)) t> « 

i? B^ 

+ ao {a\> {a\> {a> a))) + — |r (a (> (a [> a)) (> (a t> a)+ 



4! V V V jjj ' 2!2! 
^2 ^ 

+ ^ a > ((a [> (a > a)) > a) + -^Bf {a\> a)\> {{a \> a) \> a) 

+ B2BI a > (((a >a)t> a)t> a). 

It is an open question whether one can further reduce the order five term using the left pre-Lie identity. 
Moreover, it is it would be interesting to find a recursion for r2'(Aa) which already incorporates the 
pre-Lie identity. 

Acknowledgements We would like to thank A. Iserles, A. Lundervold, H. Munthe-Kaas and A. 
Munthe-Kaas-Zanna for helpful discussions and useful remarks. 
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